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1 INTRODUCTION 



The discovery of Fe-K line emission in galaxy clusters 
(Mitchell et al. 1976; Serlemitsos et al. 1977) showed that 
their X-ray emitting intracluster medium (ICM) contains 
significant amounts of processed elements created by su- 
pernovae. Different types of supernovae synthesize differ- 
ent ratios of elements. In particular, type la supernovae 
(SN la) produce large amounts of Fe and Ni. Meanwhile, 
core-collapse supernovae (SNcc) produce lighter elements 
such as O, Ne and Mg. Si-group elements (Si, S, Ar, and 
Ca) are produced by both supernova types. By measuring 
specific elemental abundances within the ICM, we can sepa- 
rate the relative contributions to the chemical enrichment 
by different types of supernovae (e.g. Werner et al. 2008; 
Bohringer & Werner 2010 and references therein). Such 
measurements place important constraints on models of 



ABSTRACT 

Using a deep (574 ks) Chandra observation of M 87, the dominant galaxy of the nearby 
Virgo Cluster, we present the best measurements to date of the radial distribution of 
metals in the central intracluster medium (ICM). Our measurements, made in 36 
independent annuli with ^250,000 counts each, extend out to a radius r ~ 40 kpc 
and show that the abundance profiles of Fe, Si, S, Ar, Ca, Ne, Mg, and Ni are all 
centrally peaked. Interestingly, the abundance profiles of Si and S - which are measured 
robustly and to high precision - are even more centrally peaked than Fe, while the 
Si/S ratio is relatively flat. These measurements challenge the standard picture of 
chemical enrichment in galaxy clusters, wherein type la supernovae (SN la) from 
an evolved stellar population are thought to dominate the central enrichment. The 
observed abundance patterns are most likely due to one or more of the following 
processes: continuing enrichment by winds of a stellar population pre-enriched by 
SNcc products; intermittent formation of massive stars in the central cooling core; 
early enrichment of the low entropy gas. We also discuss other processes that might 
have contributed to the observed radial profiles, such as a stellar initial mass function 
that changes with radius; changes in the pre-enrichment of core-collapse supernova 
progenitors; and a diversity in the elemental yields of SN la. Although systematic 
uncertainties prevent us from measuring the O abundance robustly, indications are 
that it is about 2 times lower than predicted by the enrichment models. 
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supernova explosions (e.g. Dupke & White 2000; de Plaa 
et al. 2007; Simionescu et al. 2009). 

Previous results for clusters and groups pointed to 
a centrally peaked Fe abundance distribution, coupled 
with a relatively flat O abundance profile, in cool core 
systems. Based on such measurements, it was proposed 
that the relative contribution of SN la to the enrich- 
ment of clusters increases towards their central regions 
(Finoguenov et al. 2000; Bohringer et al. 2001; Tamura 
et al. 2001; Finoguenov et al. 2002; Matsushita et al. 2003). 
The central Fe abundance peaks in cool core clusters were 
expected to form primarily by SN la with long delay times, 
as well as stellar mass loss in the cD galaxies (Bohringer 
et al. 2004; De Grandi et al. 2004). SNcc products such as 
O, Ne, and Mg were, on the other hand, produced early 
in the formation history of clusters at z ~ 2-3 and were 
thought to be well mixed throughout the ICM. Because 
SNcc also produce large amounts of Si and S, this scenario 
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predicts shallower central gradients for the distributions of 
Si and S than for Fe, resulting in Si/Fe and S/Fe ratios that 
increase with radius. 

The predicted Si/Fe gradient has, however, not been 
seen and the observed radial profiles of the Si/Fe ratios are in 
most clusters consistent with being flat, indicating that Si is 
just as peaked as Fe (e.g. Bohringer et al. 200f ; Finoguenov 
et al. 2002; Tamura et al. 2004; Sanders et al. 2004; Durret 
et al. 2005; de Plaa et al. 2006; Werner et al. 2006a; Sato 
et al. 2007, 2008; Matsushita et al. 2007). Moreover, con- 
trary to previous claims, Simionescu et al. (2009) show that 
in Hydra A and in a sample of other nearby clusters of 
galaxies the O abundance is also centrally peaked. The ob- 
served centrally peaked distribution of SNcc products sug- 
gests that the metallicity gradients form early in the history 
of clusters and persist for a long time, or that the enrich- 
ment by the winds of the evolved stellar population is signif- 
icantly more efficient than previously thought (Simionescu 
et al. 2009). Centrally peaked SNcc products are also ob- 
served in the core of the Centaurus Cluster and have been 
interpreted as being either due to continuous or intermit- 
tent star formation over the past ^8 Gyr, or due to the 
early enrichment during the formation of the central galaxy 
(see Sanders & Fabian 2006). 

This is the third in a series of papers analysing a 
deep (574 ks) Chandra observation of M87, the dominant 
galaxy of the Virgo Cluster. The first two papers (Million 
et al. 2010, hereafter Paper I; Werner et al. 2010, hereafter 
Paper II) focus on the effect of the AGN on the surround- 
ing hot plasma. This paper focuses on the history of chem- 
ical enrichment of the central regions of the Virgo Cluster. 
Many previous studies of this subject utilize XMM-Newton 
or Suzaku because of their better spectral resolution. How- 
ever, the far superior spatial resolution of Chandra, which 
allows us to better separate out obvious complications due 
to substructure (as well as the sheer number of photons al- 
ready collected with Chandra) allow us to probe the detailed 
distribution of metals with greatly improved accuracy over 
previous work. 

The structure of this paper is as follows. Sect. 2 de- 
scribes the data reduction and spatially resolved spec- 
troscopy. Sect. 3 presents our detailed abundance and abun- 
dance ratio profiles. Sect. 4 describes the implications of our 
measured radial abundance profiles on the history of chemi- 
cal enrichment through supernovae. Sect. 5 summarizes the 
results and the conclusions. 

Throughout this paper, we assume that the cluster lies 
at a distance of 16.1 Mpc (Tonry et al. 2001), for which the 
linear scale is 0.078 kpc per arcsec. 



2 DATA REDUCTION AND ANALYSIS 
2.1 Data reduction 

A detailed description of the data reduction and analysis is 
described in Paper I. The Chandra observations were carried 
out using the Advanced CCD Imaging Spectrometer (ACIS) 
between July 2002 and November 2005. The net exposure 
time after cleaning is 574 ks. Here, we use an updated ver- 
sion of the CIAO (version 4.3) software package, including 
the appropriate gain maps and updated calibration products 
(CALDB version 4.4.1). 
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Figure 1. Example spectrum and residuals from the radial anal- 
ysis. The best-fit APEC model is overplotted as a solid line and 
is a good fit to the data. The line contributions of O, Ne, Mg, 
Si, S, Ar, Ca, Fe, and Ni are labeled on the plot and are easily 
observable in spectra with ~250,000 net counts. The lower panel 
shows the residuals to the fit. 

2.2 Spectral analysis 

In order to minimize systematic uncertainties associated 
with the multi-temperature structure of the X-ray emitting 
gas, we conservatively exclude the X-ray bright arms and the 
innermost multiphase core (more details can be found in Pa- 
per I). As discussed in Paper I, beyond the innermost core 
and the X-ray bright arms, the gas can be well described by 
an isothermal model at each radius. We use partial annuli 
that vary in width from 10 to 30 arcsec, with ~250,000 net 
counts in each of a total of 36 regions. 

Background spectra were determined using the blank- 
sky fields available from the Chandra X-ray Center (see Pa- 
per I for details). Due to the high X-ray brightness of the 
target, uncertainties in the background modeling have lit- 
tle impact on the determination of the quantities presented 
here. 

The spectra have been analyzed using XSPEC (version 
12.5; Arnaud 1996). Each annulus is fit with a photoelec- 
trically absorbed (Balucinska-Church & McCammon 1992) 
single temperature APEC thermal plasma model (Smith 
et al. 2001; AtomDB v2.0.1 was used). We fix the Galactic 
absorption to 1.93 x 10 20 cm -2 (determined from the Lei- 
den/Argentine/Bonn radio H I survey; Kalberla et al. 2005). 
In order to investigate modeling uncertainties in the deter- 
mination of chemical abundances, we also repeated the spec- 
tral fits using the MEKAL thermal plasma model (Kaastra 
& Mewe 1993; Liedahl et al. 1995). All spectral fits were 
carried out in the 0.6 — 7.0 keV energy band. The extended 
C-statistic available in XSPEC, which allows for background 
subtraction, was used for all spectral fitting. 

The temperature, the normalization, and the abun- 
dances of O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni are free 
parameters for every annulus. Fig. [Tjshows an example spec- 
trum containing ~250,000 net counts with the metal line 
emission labeled. The spectra allow us to determine the Fe 
abundance to within < 3 per cent, the abundances of Si and 
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Table 1. Summary of the results of the measured radial profiles 
of the abundance ratios in the ICM. The data were fit in the 4— 
40 kpc radial range to a linear model of the form Z = a+br, where 
Z is the abundance ratio and r is the radius in kpc. Columns 
note the specific abundance ratio, best fit parameters a (in Solar 
units) and b (in units of 10 -3 Solar kpc -1 ), and the X 2 /v. fn 
cases where the APEC and MEKAL plasma codes disagree, we 
have included the best fit parameters for each code. The errors on 
parameter a are at the 68 per cent confidence level. The errors on 
the slope b are at the 68 per cent and at the 95 per cent (in the 
parenthesis) confidence level. The 95 per cent confidence upper 
and lower limits of the slope b are overplotted on each abundance 
ratio profile (Figs. EJi-f , Et-e, HJi-f). We note that only the Ni/Fe 
abundance ratio profile is not fit well by a simple linear model. 



Ratio 


a 


b 




Si/Fe 


1.316 ±0.013 


-7.5 ±0.7 (1.3) 


18.4/34 


S/Fe 


1.271 ±0.019 


-5.6 ± 1.0 (2.0) 


39.7/34 


Si/S 


1.026 ±0.019 


-1.5 ± 1.0 (2.0) 


43.4/34 


Ar/Fe 


1.00 ±0.05 


-3.1 ± 2.3 (4.7) 


34.3/34 


Ca/Fe 


1.27 ±0.06 


-1.1 ± 2.8 (5.5) 


40.1/34 


Ar/Ca 


0.78 ±0.05 


-3.7 ±2.4 (4.8) 


36.2/34 


Ne/Fe A PEC 


2.02 ±0.07 


-15.9 ± 3.3 (6.7) 


42.3/34 


Ne/Fe M EKAL 


0.80 ±0.04 


-4.4 ± 1.9 (3.8) 


43.9/34 


Mg/Fe 


0.92 ±0.02 


+3.7 ± 1.3 (2.5) 


26.5/34 


Ni/Fe APEC 


2.12 ±0.05 


-12.6 ± 2.6(5.2) 


65.9/34 


Ni/Fe MEK AL 


2.13 ±0.05 


-19.0 ± 2.5 (5.1) 


56.9/34 



S to within < 5 per cent, the abundances of Ne, Mg, and 
Ni to within < 10 per cent, and the abundances of Ar and 
Ca to within < 20 per cent statistical precision. The errors 
on the abundance profiles were determined from a Markov 
Chain Monte Carlo (MCMC) analysis. Measurement errors 
were determined from the 68 per cent confidence posterior 
distribution of the MCMC analysis. Chain lengths were at 
least 10 4 samples after correcting for burn-in. 

Abundances in the paper are given with respect to the 
'proto-Solar' values of Lodders (2003). 

2.2.1 The AtomDB v2.0.1 atomic database 

Our modelling makes use of the recently updated AtomDB 
v2.0.1 atomic database used by the APEC thermal plasma 
model implemented in XSPEC. This represents a major 
update from the previous AtomDB vl.3.2 with nearly all 
atomic data replaced^] The update of the atomic database 
affects most significantly the Fe abundance, which is on av- 
erage lower by ~ 20 per cent in v2.0.1 compared to vl.3.2. 
This change has a slight dependence upon the temperature 
of the plasma. The abundances of Si, S, Ar, Ca, Ne, Mg, 
and Ni are smaller by less than ~10 per cent in the updated 
version. The abundance ratios with respect to Fe are signifi- 
cantly higher as a result. We note that our main conclusions 
based on the centrally peaked abundances and abundance 
ratios are not sensitive to our choice of plasma code. 



3 METALLICITY OF THE X-RAY GAS 

We determine the abundance profiles for Fe, Si, S, Ar, Ca, 
Ne, Mg, and Ni and find that all elements have a centrally 
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peaked distribution. We also measure the abundance ratios 
of the individual elements with respect to Fe. In the 4-40 kpc 
radial range, we fit a linear model of the form Z = a + br 
to these measurements, where Z is the abundance ratio, r is 
the radius in kpc, a is the normalization of the linear trend, 
and b is the slope. Table [1] summarizes the best-fit linear 
relations for each of the abundance ratio profiles. For the 
Ne/Fe and Ni/Fe ratios, which have significant systematic 
uncertainties, we present the best-fit parameters for profiles 
determined using both the APEC and MEKAL plasma codes. 

The absolute abundances determined by the MEKAL 
plasma code are larger by 1-2 per cent for Si, by ~3 per 
cent for S, by ~15 per cent for Ar and Ca, and by ~20 per 
cent for Fe than the values obtained using APEC. The slopes 
of the abundance ratio profiles for these elements are, how- 
ever, consistent with those determined by the APEC code. 
The abundances of Ne, Ni, and Mg are significantly different 
when determined by the MEKAL plasma code. A discussion 
of modeling biases for the abundances of these elements can 
be found in Sect. I3.3TT1 

For gas with kT > 2.0 keV, the O abundance is ex- 
tremely difficult to determine with Chandra. At the Oviii 
line energy of 0.65 keV the effective area of the ACIS de- 
tectors is significantly affected by buildup of contamination. 
Additionally, in the lower surface brightness areas at larger 
radii the O abundance measurements are sensitive to the 
assumed Galactic foreground model. Because our O abun- 
dance measurements have large systematic uncertainties and 
may be strongly biased, we do not report their best fit val- 
ues. 

We have examined the abundance ratios separately to 
the north and south of M87. Although the overall abun- 
dances are larger to the south (see Paper I), the abundance 
ratios determined from the northern and southern sectors 
agree well with each other and the azimuthally-averaged 
analysis presented here. 

3.1 Fe, Si, and S abundances 

The top row of Fig. [2] shows the Fe, Si, and S abundance 
profiles, respectively. The bottom row shows the abundance 
ratio profiles of Si/Fe, S/Fe, and Si/S. We emphasize that 
these three elements have the most robustly determined 
abundances, with statistical uncertainties of less than 5 per 
cent. 

The Fe abundance profile (Fig. [2^) peaks at Zf c > 1.2 
Solar in the central regions. The Si and S abundance profiles 
(Fig. [2}>c) peak at a slightly larger central value of Zsi ; s ~ 
1.5 Solar. These profiles then exhibit steady declines to ~ 0.6 
Solar by r ~ 25 kpc. 

A significant enhancement, or 'bump', in the Fe and Si 
(and possibly S) abundances is seen at r ~30 kpc. This is 
approximately the radius at which the bright X-ray arms 
terminate. As discussed in Paper I, this bump may be due 
to the uplift of cool, metal-rich material in the wake of buoy- 
antly rising radio bubbles. 

Most importantly, we observe, for the first time, a radi- 
ally decreasing trend in the Si/Fe and S/Fe abundance ratios 
(Fig.[5Jl-e). They both peak at ~ 1.3 Solar near the core and 
decline to ~ 1.0 Solar by r ~ 35 kpc. Both the Si/Fe and 
S/Fe profiles are well described by a linear decline as a func- 
tion of radius. As seen in Table [TJ both profiles have similar 
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Figure 2. Top row: Abundance profiles of Fe (left; see also Paper I), Si (middle), and S (right). Bottom row: Abundance ratio profiles 
of Si/Fe (left), S/Fe (middle), and Si/S (right). Fe, Si, and S have the best determined abundances. The 95 per cent lower and upper 
limits of the slopes of the radial distributions of the abundance ratios are overplotted. 



slopes to one another and are significantly steeper than a 
constant ratio (at the 5-1 la level). The Si/Fe and S/Fe ra- 
tios of ~0.9 Solar measured at radii between 50-200 kpc 
with Suzaku (Simionescu et al. 2010) indicates that beyond 
r ~ 35 kpc these abundance ratio profiles flatten. 

The Si/S abundance ratio profile (Fig. [2J) is an interest- 
ing cross check upon the determination of these abundances. 
Because both elements are created in roughly equal quan- 
tities by both SNcc and SN la (see Iwamoto et al. 1999; 
Nomoto et al. 2006), the Si/S abundance ratio is insensitive 
to the relative fraction of supernovae that explode as SNcc 
and SN la. Any trend in the Si/S abundance ratio profile is 
instead primarily affected by a change in the average yields 
of SNcc and SN la as a function of radius. The Si/S abun- 
dance ratio is consistent with being flat as a function of ra- 
dius and its observed mean value is Zsi/Zs = 1.002 ± 0.009 
Solar. 

The systematic uncertainties involved in the modelling 
of Si and S lines are relatively small due to the favorable 
location of their lines in the 2-3 keV range, where no resid- 
ual Fe-L emission is present. Possible biases due to multi- 
temperature structure are minimized by our choice of ex- 
traction regions which avoid clear density and tempera- 
ture inhomogeneities and have an approximately isothermal 
azimuthal temperature structure (see Paper I). Nonethe- 
less, we performed detailed spectral simulations of multi- 
temperature plasmas to investigate possible biases due to 
residual structure and projection effects. We mixed APEC 



thermal models with slightly different temperatures and 
metallicities, where the cooler plasma is more enriched than 
the hotter one. The simulated gas mixtures span a grid 
of possible projection effects and unresolved temperature 
structure. We fitted these simulated spectra with a single 
temperature model. The potential bias in Si/Fe, S/Fe and 
Si/S ratios was always less than 10 per cent. 

Furthermore, we note that systematic uncertainties re- 
lated to the effective area calibration of the ACIS detectors 
around the Si edge would be likely to affect the Si and S 
abundance profiles differently. Therefore, the striking simi- 
larity of the Si/Fe and S/Fe profiles strongly indicates that 
the detection of these trends is robust. 

3.2 Ar and Ca abundances 

The top row of Fig. [3] shows the Ar (left) and Ca (right) 
abundance profiles. The bottom row shows the radial dis- 
tributions of the Ar/Fe (left), Ca/Fe (middle), and Ar/Ca 
(right) abundance ratios. Because the Ar and Ca lines have 
small equivalent widths, their abundances are more difficult 
to measure precisely than Fe, Si, and S: our data enable 
us to measure the Ar and Ca abundances to within 20 per 
cent statistical uncertainty. Our spectral fits reveal that the 
Ar and Ca abundances are relatively insensitive to small 
changes in the temperature and Fe abundance. Therefore, 
systematic uncertainties on these measurements are likely 
small. 
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Figure 3. Top row: abundance profiles of Ar (left) and Ca (right). Bottom row: abundance ratio profiles of Ar/Fe (left), Ca/Fe (middle), 
and Ar /Ca (right) . The 95 per cent lower and upper limits of the slopes of the radial distributions of the abundance ratios are overplotted. 



The Ar and Ca abundance profiles (Fig. "Ji-b) are 
also centrally peaked. The Ar and Ca abundances peak at 
Z\ T ~ 1.2 and Zc a ~ 1.4 Solar, respectively. They decrease 
to Zat ~ 0.5 and Zc& ~ 0.6 Solar respectively by r ~ 35 kpc. 
As with Fe and Si, both profiles show a marginal, but plausi- 
ble increase in abundance at r ~ 30 kpc that may be due to 
the uplift of cool, metal-rich material, as described in Paper 
I. 

Both the Ar/Fe and Ca/Fe abundance ratio profiles 
(Fig.|3J>d) are consistent with constant values of Zay/Zys = 
0.94 ± 0.02 and Z Cll /Z Fc = 1.25 ± 0.03 Solar, respectively. 

Similarly to the Si/S abundance ratio profile, the Ar/Ca 
abundance ratio profile (Fig. is also an interesting cross- 
check upon the determination of these abundances. Like Si 
and S, Ar and Ca are also created in similar quantities by 
SNcc and SN la. Therefore, the Ar/Ca abundance ratio pro- 
file is insensitive to the relative number of supernovae that 
explode as SNcc and SN la. Instead, the Ar/Ca abundance 
ratio is primarily sensitive to changes in the average yields 
of SN la. The Ar/Ca abundance ratio is consistent with be- 
ing flat as a function of radius and its observed mean value 
is Zat/Zcz = 0.71 ± 0.02 Solar. 

3.3 Ne, Mg, and Ni abundances 

The top row of Fig. ["] shows the abundance profiles of Ne 
(left), Mg (middle), and Ni (right), respectively. The bot- 
tom row shows the abundance ratio profiles of Ne/Fe (left), 
Mg/Fe (middle), and Ni/Fe (right). There are significant 



systematic uncertainties present in the determination of the 
Ne, Mg, and Ni abundances (see Sect. 13.3711 . 

The inferred abundance profiles of Ne, Mg, and Ni are 
all centrally peaked at Zn c ~ 2.3 Solar, Z-Mg, ~ 1.2 Solar, 
and ^Ni ~ 2.8 Solar, respectively. The profiles decline to 
Znc ~ 1-0 Solar, Zm s ~ 0.6 Solar, and Zm ~ 1.0 Solar 
by r ~ 25 kpc. Although there are differences between the 
abundances measured with APEC and MEKAL, strong cen- 
tral peaks in the abundance profiles of Ne, Mg, and Ni are 
found with both codes. These profiles also show the enhance- 
ments at ~30 kpc. The Ne/Fe abundance ratio determined 
by both plasma codes is centrally peaked. 

The Mg/Fe ratio is marginally consistent with being flat 
as a function of radius and its mean values is Zyi^/Z-pe = 
0.988 ± 0.012 Solar. This is in agreement with previous re- 
sults once the differences between plasma codes are taken 
into account, (see e.g. Matsushita et al. 2003, Simionescu 
et al. 2010). 

The Ni/Fe results derived using both plasma codes are 
poorly fit to a linear model. The reduced x 2 s are high, 
X 2 /v ~ 60/34. Primarily this is due to increased scatter 
in the Ni measurements at large radii (r > 30 kpc). Within 
the radial range of 5 < r < 30 kpc and using the APEC 
code the Ni/Fe abundance ratio is fit well by a single, con- 
stant value of Zn\/Zyc = 1.92 ±0.02 Solar. However, a large 
drop in the Ni/Fe abundance ratio at radii r > 30 kpc is 
reproduced by both plasma codes. 
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Figure 4. Top row: abundance profiles of Ne (left), Mg (middle), and Ni (right). Bottom row: abundance ratio profiles of Ne/Fe (left), 
Mg/Fe (middle), and Ni/Fe (right). The Ne, Mg, and Ni lines are blended with the Fe-L line emission. Therefore, significant modeling 
bias may be present in the determination of the abundances of these elements. We note, however, that the abundances of Ne, Mg, and 
Ni are all centrally peaked, regardless of the choice of plasma code. The 95 per cent lower and upper limits of the slopes of the radial 
distributions of the abundance ratios are overplotted. 



3.3.1 Modelling uncertainties in Ne, Mg, and Ni 
measurements 

The spectral resolution of the CCD type detectors on Chan- 
dra does not allow us to resolve the individual Ne and Ni 
lines within the Fe-L complex. The Mg lines are also blended 
with Fe-L line emission. Therefore, the Ne, Mg, and Ni abun- 
dances dependent critically on the modeling of the Fe-L com- 
plex. 

The Ne, Mg, and Ni abundances determined by the 
MEKAL and APEC plasma codes disagree, with the MEKAL 
Mg abundances being a factor of two lower than the values 
determined with APEC. Both plasma models remain incom- 
plete and therefore the results must be treated with caution. 

Previous results from the Reflective Grating Spectrom- 
eter (RGS) aboard XMM-Newton, which can resolve sepa- 
rately the Ne lines within the Fe-L complex, reveal that the 
Ne/Fe abundance ratio is above the Solar value (see Werner 
et al. 2006b), consistent with our Ne abundance measure- 
ment. 

Incorrect ionization balance for Ni in the plasma models 
can lead to temperature dependent biases in its abundance 
determination (see Bohringer & Werner 2010). At kT ~ 2.1 
keV, the Ni abundances determined with APEC are 20 per 
cent lower than those determined using MEKAL. This dif- 
ference decreases linearly with temperature. At kT ~ 2.7 



keV, the best fit Ni abundances determined using the two 
plasma codes are consistent. Our Ni abundances are deter- 
mined from their L-shell lines. Using additional spectral data 
between 7.0 — 8.0 keV to include the K-shell lines does not 
improve our determination of the Ni abundance. 



4 IMPLICATIONS FOR SUPERNOVA 
ENRICHMENT MODELS 

The elements observed in the X-ray emitting gas of clusters 
of galaxies reveal the integrated chemical enrichment history 
by supernovae. The standard picture of chemical enrichment 
in clusters of galaxies for at least the past 12 years has been 
an early enrichment by SNcc products that are now well 
mixed in the ICM, in combination with a later contribution 
by SN la, which have longer delay times and primarily en- 
rich the region surrounding the cD galaxy. Our data clearly 
disagree with this picture. Both Si and S show central abun- 
dance peaks that are larger than that of Fe (see Figs. d£-f, 
[4}i). The abundances of other elements (Ar, Ca, Ne, and Mg) 
show central abundance peaks as well (see Figs. [3ji-d, |4^). 
These results force us to rethink our models of the chemical 
enrichment in clusters of galaxies. 

Here we discuss possible chemical enrichment scenar- 
ios that may explain the observed abundance ratio profiles. 
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These include centrally peaked SNcc products due to stel- 
lar winds, intermittent star formation, and early enrichment 
of the low entropy gas. We also discuss how radial changes 
in the stellar initial mass function, and pre-enrichment of 
SNcc progenitors, or a diversity in the population of SN la 
can affect the observed distributions of metals. 



1.16 x 10 ergs, which was shown to fit best the properties 
of the Tycho supernova remnant (Badenes et al. 2006). 

The predicted Mg/Fe abundance ratio profile disagrees 
with the measurements. The Mg abundance measurements 
are, however, affected by systematic uncertainties arising 
from the incomplete modeling of the Fe-L shell transitions. 



4.1 Centrally peaked SNcc products 

The most straight forward interpretation of the centrally 
peaked Si/Fe and S/Fe abundance ratio profiles, is that there 
is a dominant contribution of SNcc products to the enrich- 
ment of the lowest entropy gas in the central regions of the 
galaxy. The contribution of SN la products rises with in- 
creasing radius, out to r ~ 35 kpc. 

We calculate the fraction of all supernovae that explode 
as SN la as a function of radius using the equation 



(Mi \ 
\M F J 



_ ^ /M Z ;SNIa + (1 - /)M Z;S ] 

Fe J obs ~ /M Fc;S NIa + (1 - /)Mp c;S 



(1) 



where (Mz/MF C ) obB is the mass ratio of element Z with 
respect to Fe converted from our abundance measure- 
ments using the proto-Solar values by Lodders (2003), / = 
A r sNia/(A^sNia + AfsN cc ) is the fraction of all supernovae that 
explode as SN la, and JWz ; sNia and Mz ; sn cc are the aver- 
age theoretical yields for SN la and SNcc, respectively. For 
SN la, we use the yields of the WDD2 delayed-detonation 
model by Iwamoto et al. (1999), except where explicitly 
stated. Because Si is the second most reliably measured 
chemical element after Fe, we calculate the number fraction 
of SN la, f(r), using the Si/Fe ratio. 

The average theoretical yields of SNcc are calculated 
using the values by Nomoto et al. (2006) under the assump- 
tion of the Salpeter initial mass function (IMF). In detail, 



Z;SN CC = 



5OM 
lOMp, 



Mz(M)M~ a dM 



J II 



M~<*dM 



(2) 



where Mz ; sn cc is the theoretical yield of element Z for the 
weighted-average of core-collapse supernovae, Mz(M) are 
the atomic yields as a function of stellar mass for core- 
collapse supernovae, and a — 2.35 for the Salpeter IMF. 

We calculate a central fraction of SN la of / = 0.10 at 
radius r ~ 50 arcsec (3.9 kpc). This fraction rises to f — 0.19 
at radius r ~ 500 arcsec (38.8 kpc). Based on these fractions, 
we can predict the radial profiles of the abundance ratios for 
other elements. The results are shown in Fig. [5] 

The S/Fe and Si/S abundance ratio profiles predicted 
by the central increase in the relative fraction of SNcc agree 
well with our observed profiles. While the observed slopes 
of the Ca/Fe and Ar/Fe profiles are also consistent with the 
centrally peaked SNcc enrichment, their normalizations are 
higher than those predicted by the model. Our Ca/Fe abun- 
dance ratio is, however, consistent with the average value 
measured in 22 clusters of galaxies by de Plaa et al. (2007) 
once systematic differences in the plasma codes used are 
taken into account. These authors showed that the measured 
abundance ratios can be well fit with a one-dimensional 
delayed-detonation SN la model, calculated on a grid in- 
troduced in Badenes et al. (2003), with a deflagration to 
detonation density of 2.2 x 10 7 g cm -3 and kinetic energy of 



4-1.1 Centrally peaked SNcc products due to intermittent 
star-formation 

If both the low entropy gas in the cluster core and the bulk of 
the surrounding ICM were enriched around the same time, 
with a similar mixture of SNcc and SN la, then the expected 
Si/Fe and S/Fe ratios would be constant with radius. How- 
ever, out to a radius of r = 35 kpc, the observed Si/Fe and 
S/Fe ratios decreases from ~1.3 to ~1.0 Solar, before stay- 
ing flat out to at least a radius of 200 kpc (Simionescu et al. 
2010). In order to produce the observed Si/Fe peak on top 
of this flat profile, an additional 1.6 x 10 7 SNcc would have 
to explode in the centre of M 87. This is about 10 per cent 
of all SNcc that are required for the chemical enrichment of 
the innermost r < 40 kpc region. Assuming a Salpeter IMF, 
such number of supernovae require ~ 10 9 Mq of star forma- 
tion. Because of the continuing production of Fe by SN la in 
the cD galaxy, this value is a lower limit. The SNcc enrich- 
ment would also have to occur on a time scale shorter than 
the production of Fe by SN la. While, there is no evidence 
for current star formation in M 87, and the 95 per cent confi- 
dence upper limit on radiative cooling from the coolest ICM 
phase is 0.06 Mq yr _1 (Werner et al. 2010), we can not rule 
out past intermittent star forming episodes in the central 
regions of the galaxy. Star formation, at rates reaching sev- 
eral tens and up to hundreds of solar masses per year, has 
been observed in the brightest cluster galaxies of some other 
cooling core clusters (e.g. McNamara et al. 2006, Ogrean et 
al. 2009, Ehlert et al. 2011). Therefore, several intermittent 
star forming episodes at ~10 Mq yr _1 lasting for a total of 
~2x 10 s yr would not be surprising. 



4- 1.2 Centrally peaked SNcc products due to stellar winds 

Stellar mass loss is an important source of metals in the hot 
gas surrounding giant elliptical galaxies like M 87. The ma- 
terial that formed the current stellar population of M 87 was 
most likely pre-enriched primarily by SNcc products. Some 
of these metals are then returned into the ICM/ISM by stel- 
lar winds. Assuming a single-age passively evolving stellar 
population with a Salpeter initial mass function, Ciotti et 
al. (1991) predict a stellar mass loss rate of 

M*(t) w 1.5 x lO -11 !®^ 1-3 



(3) 



where Lb is the present-day B-band luminosity in units of 
Lbq (L b = 8.4 x 10 10 £ B q for M 87; Gavazzi et al. 2005) 
and tis is the age of the stellar population in units of 15 Gyr 
(the formula is valid in the range from ~0.5 to over 15 
Gyr). Integrating this equation under the assumption that 
the current stellar population is 10.5 Gyr old (formed at 
z — 2.1), we obtain a gas mass contribution from stellar 
winds of 1 x 10 11 Mq during the past 10 Gyr. Assum- 
ing that the material from which the current stellar pop- 
ulation of M 87 formed was pre-enriched to conservatively 
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Figure 5. Predicted abundance ratio profiles using the calculated fraction of SN la, f(r), from the Si/Fe abundance ratio profile with the 
measured data points at 4 kpc and 35 kpc over-plotted. The predicted abundance ratios of S/Fe and Si/S agree well with observations. 
However, the predicted radial distributions of the abundance ratios of the other elements do not match the observations well. 



Zsi ~ 0.5 Solar, with a Si/Fe abundance ratio of 1.5 Solar, 
the total mass of Si returned to the ISM/ICM by the stel- 
lar winds is 4.3 x 10 7 M®. The total mass of Fe returned is 
4.8 x 10 7 Mq. Under these assumptions, the total mass of Si 
produced by stellar winds in excess of a Si/Fe ratio of 1 Solar 
is 1.4 x 10 7 M Q . The observed total Si mass in excess to a 
flat profile with Si/Fe=l Solar around M 87 is 3.6 x 10 6 M@. 
Therefore, despite all the uncertainties in the estimates of 
the metal mass loss, the excess Si observed around M 87 
could most likely be produced by stellar winds. Because the 
initial starbursts, that enriched the material from which the 
current stellar population formed, produced predominantly 
SNcc products, and the Mg2 index indicates that the stel- 
lar population of M 87 is enriched to more that 1 Solar 
(Kobayashi & Arimoto 1999, Matsushita et al. 2003) this 
scenario can plausibly produce centrally peaked abundance 
profiles of SNcc products. 

4-1.3 Centrally peaked SNcc products due to strong early 
enrichment of the low entropy gas 

Another possible mechanism contributing to the observed 
centrally peaked distribution of SNcc products is strong 
early enrichment of the lowest entropy X-ray emitting gas 
and inefficient mixing of this material with the surround- 
ing ICM. If the lowest entropy X-ray emitting gas currently 
seen in the Virgo Cluster core was at high redshift located in 
the environments of massive galaxies during their epoch of 
maximum star formation, then this material may have be- 
come enriched significantly by SNcc products. As the cluster 
formed, this low entropy SNcc enriched gas will naturally 
sink and accumulate at the base of the cluster potential. 
Assuming that it does not become well mixed with the sur- 
rounding ICM as it does so and it does not cool out of the 
hot phase, this may lead to the observed peak in metal abun- 
dances. 

4.2 Other mechanisms affecting the radial profiles 
of abundance ratios 

We explore other possible explanations for the observed 
abundance ratio profiles. Each of these mechanisms predicts 
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Figure 6. Predicted abundance ratios as a function of the slope 
of the stellar initial mass function. The value 2.35 corresponds 
to the Salpeter value. Flatter mass functions, i.e. larger relative 
fractions of massive stars exploding as SNcCi predict a rise in 
Si/Fe and S/Fe ratios as we observe in the centre of the galaxy. 
Based on SNcc yields by Nomoto et al. (2006) and assuming that 
SN la (WDD2 model by Iwamoto et al. 1999) make up 15 per cent 
of all supernovae. 

a central enhancement in the Si/Fe, S/Fe, Ar/Fe, and Ca/Fe 
ratios and each of them might, to some degree, contribute 
to the observed radial trends. 

4-2.1 Changing IMF as a function of radius 

A stellar IMF, which is flatter (has a smaller a in Equation 
2, therefore producing more massive stars) in the central 
regions of the cluster would produce a central increase of Si- 
group elements. The ratios of chemical elements produced by 
SNcc are a strong function of the mass of the progenitor. To 
examine the effect of the IMF on the predicted abundance 
ratios, we vary its slope a. 

Fig. [S] shows the predicted abundance ratios assuming 
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Figure 7. Predicted abundance ratios as a function of the initial 
metallicity of SNcc progenitors. The absolute metallicity of Zi = 
0.02 represents the Solar metallicity. Based on SNcc yields by 
Nomoto et al. (2006) and assuming that SN la (WDD2 model by 
Iwamoto et al. 1999) make up 15 per cent of all supernovae. 
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Figure 8. Predicted abundance ratios for a variety of SN la 
yield models. We considered the WDD1, WDD2, WDD3, and 
W7 models of Iwamoto et al. (1999). For SNcc we use the yields 
by Nomoto et al. (2006) and we assume that SN la make up 15 
per cent of all supernovae. 



a steepening IMF as a function of radius and a radially con- 
stant relative fraction of SN la, / = 0.15. This fraction was 
chosen to lie within the range suggested by the Si/Fe abun- 
dance ratio profile and it only affects the normalization of 
the predicted profiles. The explanation of the observed range 
of values in the radial profiles of the Si/Fe and S/Fe ratios 
by a radial trend in the IMF would require extremely steep 
IMF at larger radii. Moreover, this scenario would produce 
very large central increases of the Ne/Fe and Mg/Fe ratios, 
which we would clearly detect even given the systematic un- 
certainties in the Ne and Mg abundance determinations. 

While we cannot rule out a small radial trend in the 
stellar IMF, it cannot be responsible for the large observed 
ranges in the abundance ratios. 

4-2.2 Radial trends in the SNcc pre- enrichment 

We examine the effects of a possible radial trend in the SNcc 
pre-enrichment, i.e. the metallicity of the stars that produce 
the supernovae on the radial profiles of abundance ratios in 
the ICM. 

Fig. [7] shows the predicted abundance ratios assuming 
that the metallicity of the SNcc progenitors is increasing 
as a function of radius. The yields of SNcc as a function 
of the initial metallicity are from Nomoto et al. (2006) and 
the absolute metallicity of Zi = 0.02 is equal to Solar. We 
assume a constant relative fraction of SN la (/ = 0.15) and 
the Salpeter IMF. The plot shows that SNcc with a lower 
initial metallicity produce higher Si/Fe and S/Fe ratios. 

Enrichment from infalling, low-entropy systems, which 
may be dominated by SNcc with relatively low metallicity 
progenitors during starbursts, could have contributed to the 
observed central peaks in the abundance ratios. 

We note that SNcc with a very low initial metallic- 
ity (Zi < 0.004) would produce large Si/Fe, S/Fe, Mg/Fe, 
Ca/Fe, Ar/Fe, Si/S, and Ar/Ca ratios and small Ne/Fe and 



Ni/Fe ratios that do not match our observations. Therefore, 
bulk of the metals observed in the ICM could not have been 
produced by extremely low metallicity stars. 

4-2.3 Diversity in the SN la population 

There is a growing evidence for a diversity in SN la ex- 
plosions (see e.g. Sullivan et al. 2006; Pritchet et al. 2008; 
Mannucci et al. 2006). While a population of brighter SN la 
with a slow luminosity decline is more common in late-type 
spiral and irregular galaxies with recent star formation (in- 
dicating a short delay time between their formation and the 
explosion), a fainter and more rapidly decaying population 
of SN la is more common in early-type galaxies (Hamuy 
et al. 1996; Ivanov et al. 2000). This diversity should be re- 
flected in the abundance yields with the brighter SN la pro- 
ducing more Ni and less Si-group elements than the fainter 
ones. 

Based on early work with XMM-Newton, Finoguenov 
et al. (2002) argue that the diversity in the SN la population 
would explain the distribution of chemical elements in the 
Virgo Cluster. 

The variation of the peak brightness, which correlates 
with the production of 58 Ni and anti-correlates with the pro- 
duction of Si-group elements, can also be explained in the 
framework of the delayed detonation models by a variation 
of the defiagration-to-detonation transition density (transi- 
tion from subsonic to supersonic flame velocities). Fig. [8] 
presents the expected yields for a variety of SN la explo- 
sion models from Iwamoto et al. (1999), with the WDD1, 
WDD2, WDD3 and W7 models on the x-axis. We assume 
a constant relative fraction of SN la, / = 0.15. The W7 
model represents a pure deflagration explosion mechanism. 
The WDD models represent delayed-detonation explosions 
and the last digit indicates the density at which the flame ve- 
locity becomes supersonic (defiagration-to-detonation tran- 
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sition density) in units of 10 7 g cm~ 3 . This transition density 
is likely dependent on the composition of the progenitor (see 
Jackson et al. 2010). 

Matching our observed Si/Fe and S/Fe profiles, given 
the existing SN la yields, would require that the centre of 
the galaxy has been enriched almost solely by WDD1 super- 
novae while the outer regions almost solely by WDD3. The 
contribution by SN la with longer delay times and larger 
Si/Fe ratio would be the largest in the centre of the galaxy, 
but any realistic enrichment scenario predicts an enrichment 
by a mixture of different types of SN la at all radii. 

This model predicts a large central increase in the Ar/Fe 
and Ca/Fe abundance ratios that seem to be in conflict with 
the observed relatively fiat profiles. The predicted Mg/Fe 
profile is relatively fiat in agreement with observations. The 
strongest prediction of this model is the ~ 30 per cent rise in 
the Ni/Fe abundance ratio with the increasing radius. Our 
observed Ni/Fe profile, which is unfortunately dominated 
by systematic uncertainties, suggests a relatively flat radial 
distribution. 

4.3 The O abundances 

The O/Fe ratio of 0.60 ±0.03 Solar determined in the centre 
of M 87 with the XMM-Newton Reflection Grating Spec- 
trometers (Werner et al. 2006b; these data resolve the O vm 
line and the individual lines of the Fe-L complex) is signifi- 
cantly lower than the values predicted by the proposed en- 
richment scenarios. These best fit O/Fe ratios are consistent 
with the values determined using the CCD type detectors 
on XMM-Newton (Matsushita et al. 2003) . Either the mea- 
surements strongly underestimate the O abundance or some 
key aspect of the chemical enrichment of the hot 1CM/ISM is 
not understood. Using the recently updated AtomDB atomic 
database, the best fitting O/Fe ratios are 50 per cent larger 
compared to the previous version, indicating that at least 
part of the discrepancy might be a modeling issue. Fur- 
thermore, all of our proposed enrichment scenarios are in- 
compatible with the rising O /Fe abundance profile reported 
from XMM-Newton (Bohringer et al. 2001, Finoguenov et 
al. 2002, Matsushita et al. 2003). However, measurements of 
the Oviii line emission at E ~ 0.65 keV with CCD detec- 
tors suffer from significant systematic uncertainties due to 
a combination of limited spectral resolution, residual gain 
uncertainties, coupled with incomplete modelling of the de- 
tector oxygen edge and possible incomplete subtraction of 
the O vm line emission from the Galactic foreground (which 
could bias the O abundance measurement in the outskirts of 
M 87 high) . These systematic uncertainties force us to treat 
all current O abundance measurements with caution. More 
robust measurements of the O/Fe profile will be possible 
with the calorimeters on the Astro-H satellite. 



5 SUMMARY 

Using a deep (574 ks) Chandra observation of M 87, we per- 
formed the best measurements to date of the radial distri- 
butions of metals in the ambient central ICM of the Virgo 
Cluster. We conclude that: 

• The abundance profiles of Fe, Si, S, Ar, Ca, Ne, Mg, 
and Ni are all centrally peaked. 



• The abundance profiles of Si and S are more centrally 
peaked than Fe, challenging the standard picture of chem- 
ical enrichment in galaxy clusters, wherein SN la products 
are thought to dominate the central enrichment. Rather, 
despite a negligible current star formation rate in M 87 and 
the continuing enrichment by SN la, the integrated relative 
contribution of core collapse supernovae to the enrichment is 
higher in the central low entropy core than in the surround- 
ing ICM. The observed abundance patterns are most likely 
due to one or more of the following processes: continuing 
enrichment by winds of a stellar population which has been 
pre-enriched mainly by SNcc products; intermittent forma- 
tion of massive stars in the central cooling core on time 
scales shorter than the continuing enrichment by SN la; or 
strong early SNcc enrichment of the lowest entropy X-ray 
emitting gas that is subsequently not well mixed and does 
not cool out of the hot ICM phase. 

• Other processes, such as stellar initial mass function 
that changes with radius, changes in the pre-enrichment of 
core-collapse supernova progenitors, and diversity in the el- 
emental yields of SN la, might have also contributed to the 
observed radial profiles. 

• Although systematic uncertainties prevent us from 
measuring the O abundance robustly, indications are that 
it is about 2 times lower than predicted by the enrichment 
models. 
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